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Towards the application of the Maximum Entropy Method to finite 
temperature Upsilon Spectroscopy * 

M. Oevers^ with C. Davies'^ and J. Shigcmitsu'^ 

'^Dept. of Physics and Astronomy, University of Glasgow G12 8QQ, U.K. 

''The Ohio State University, Columbus, Ohio 43210 

According to the Narnhofer Thirring Theorem fll interacting systems at finite temperature cannot be described 
by particles with a sharp dispersion law. It is therefore mandatory to develop new methods to extract particle 
masses at finite temperature. The Maximum Entropy method offers a path to obtain the spectral function of a 
particle correlation function directly. We have implemented the method and tested it with zero temperature T 
correlation functions obtained from an NRQCD simulation. Results for different smearing functions are discussed. 
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1. Introduction 

With J/^-suppression being one of the key 
probes for the quark gluon plasma |2| a nonper- 
turbative understanding of heavy quarkonia at 
temperatures above the transition temperature 
is important. The extraction of particle masses 
at finite temperature is aggravated on a practi- 
cal level by the compactified Euclidean time di- 
rection and on a more fundamental level by the 
Narnhofer Thirring theorem. One has to employ 
methods that do not make any assumptions about 
the spectral structure. A zero momentum Eu- 
clidean current current correlator G{t) has the 
following spectral representation M: 



G{t) = - 



cosh((§-r)u;) 
sinh('|a') 



p{uj)duj 



(1) 



The (temperature dependent) spectral function 
p{io) contains all the real-time physics informa- 
tion we are interested in. Indeed the Fourier 
transform of the retarded correlator is given in 
terms of the spectral function as 



Gr{l^) 



1 



piu') 



n J _nn ijJ + If- — IjJ 



-duj' 



Extraction of p{uj) from Eq.dy) by inversion is nu- 
merically an ill-posed problem Gj and the Maxi- 

*Work supported by the European Community TMR Pro- 
gram ERBFMRX-CT97-0122 



mum Entropy method can be seen as a regular- 
ization of this ill-posed problem, but has in fact 
a deeper justification from Bayesian statistics. 

2. The Maximum Entropy Method 

The Maximum Entropy method (MEM) is a 
well known technique for image reconstruction 
and has been successfully applied in astronomy 
and condensed matter physics. For a review see 
0. MEM is based on Bayesian methods of in- 
ference which in turn centers around Bayes theo- 
rem of conditional probabilities, which in our case 
reads: 



P[p\G,I]^P[G\p,I]P[p\I] 



(2) 



where p is the spectral function, G are the data 
for the correlator and I is any a priori information 
that is relevant to the problem. P [G\p, I] is called 
the likelihood and is proportional to exp(— L) for 
a large number of measurements with 



-k 



a 



= i;T.'^F{n)-G{r,)GrJ{F{r, 



G{t,)) 



jj being the covariance matrix and F{t) — 
J K{T,uj)p{uj)duj is the 'fit function' in terms of 
the spectral density p{uj) and the kernel K{t,uj) 
defined by Eq(|l|). P [/9|/] is the prior probability 
for the spectral function. Using Bayesian lines of 



thought I2I one can show that the prior probabil- 
ity is of entropic form exp(— aS*). 



[p{uj) - miuj) - p{uj) log{ 



m(Lu) 



)]du 



m{uj) is a default model with respect to which we 
measure the entropy of the spectral function. The 
freedom to choose a default model can be used 
to incorporate further knowledge, but trustable 
results should not depend on miuj). We choose 



m{ijj) 



which is the perturbative high en- 



ergy form of mesonic spectral functions 0j8| . The 
real and positive parameter a controls the rela- 
tive weight between the entropy and the the likeli- 
hood. In the algorithm used here to determine the 
spectral function ||9|] a is eliminated by marginal- 
ization. A detailed analysis shows that one can 
determine the probability P[a\G] from the data. 
And the final result is then given as a weighted 
average over Pa{uj) which is the spectral function 
that maximises Q — aS — L, i.e. maximises the 
conditional probability P [p|G,/] in Eq.(||): 

p{lo) = I pa{ijj)P[a-\G\da 

The maximization of Q over the space of p makes 
use of a Singular Value Decomposition of the ker- 
nel K{t,u!), by expressing p in terms of the sin- 
gular vectors of K. In this way the algorithm 
chooses the appropriate degrees of freedom of 
which there are, by construction, always fewer 
than the number of timeslices. MEM therefore 
makes no assumption about the spectral shape 
except for that which is dictated by the the dis- 
cretisation of the kernel. 

3. Testing MEM with T data 

Since we are ultimately interested in studying 
the melting of T and J/^ above the critical tem- 
perature, we have tested the method with data 
from a precision T spectroscopy study at zero 
temperature using NRQCD reported in [|l^. We 
have analyzed local and smeared correlation func- 
tions in the ^5i-channcl with Richardson poten- 
tial radial wave functions for the smearings. The 
smearings at source and sink have to be identical 
for Eq.p) to hold. Fig.n^ shows our results. 
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Figure 1. Spectral functions p{ijj)/uo'^ for (top to 
bottom) the local-local correlator, ground-state 
smearing, 1st excited state smearing, 2nd excited 
state smearing. The vertical lines indicate the 
error-bands of the results of lloi . 



In all graphs we plot /(w) = p{lo)/uP to divide 
out the assumed default model m{ijj) = O.lw^. 
The integration in Eq.(|l|) is performed up to 
^max = 7.5 > 27r. The results do not depend on 
i^max as long as cumax > Stt, the maximum lattice 
momentum of the meson. We only show the spec- 
tral function up to w ~ 1.2, since /(w) is virtually 
zero beyond this value except for w close to tUmax 
when /(w) tends to 0.1, the default model value. 
The data do not contain enough information to 
constrain the spectral function at such high mo- 
menta. The results also do not depend on the 
discretisation dui as long as duj is smaller than 
the finest structure that is contained in the data. 
These results have been produced with du = 0.01. 
We have also checked the independence of our re- 
sults from the default model. We have varied the 
prefactor of w" between 0.002 and 5.0 and the 
power n of cu" between and 3. 

3.1. Discussion of the results 

The local-local correlator shows a clear ground- 
state peak at the location expected from the anal- 
ysis in |10(| . There is some indication for spectral 
strength of excited states, but the shape is not 
very pronounced. Similar behaviour has been re- 
ported in |0,|[ . What is new here is the analysis of 
smeared correlators. Although a smeared correla- 
tion function probes the same quantum numbers 
as the local correlator, the shape of the spectral 
function of a smeared operator will in general be 
different from the local spectral function. How- 
ever the position of a peak in the spectral func- 
tion should not be affected by smearing. The sec- 
ond plot of Fig.|l| shows the result for ground-state 
smearing and indeed one finds a clear peak with 
no indication of any excited states. The next plot 
is the smearing for the first excited state. The 
peak position is consistent with the results from 
the standard analysis. There is no contamination 
from the ground-state and only a slight indica- 
tion of higher excited states. For the last plot a 
comment is in order. We have not said anything 
so far about the statistical significance of our re- 
sults and if we take the last plot at face value it 
presents us with a problem. The lower of the two 
peaks lies lower than the result for the first ex- 
cited state. Either the 1st excited state smearing 



did in fact project onto a higher state and this is 
the true first excited state or something is prob- 
lematic with the 2nd excited state smearing. The 
second peak which we expect to be the second 
excited state lies higher than the result from the 
correlated matrix fit of |l^ . The large error band 
from the standard analysis already indicates that 
this state is difficult to analyse. The error analysis 
described in the next section indicates that both 
peak positions are in fact compatible with the re- 
sults of ^^ and that the data for this smearing 
are probably not good enough to pin down this 
state. One expects that higher excited states be- 
come more and more visible as one increases Nr, 
the number of points in the time direction. At 
finite temperature it will therefore be advisable 
to use anisotropic lattices in order to have large 
Nr > 24 at still manageable spatial lattice sizes. 

3.2. Rotated smearings 

Suppose one has a set of smearings 
|0} , |1} , |2) .... Measuring all cross correlators 
(0|0) , (0|1) , (1|0) ..., one can construct all corre- 
lators with smearings that can be expressed as 
linear superpositions of the original smearings: 
1^) = ao |0) + ai |1) + 02 |2} + .... The cor- 
relator G{t) = (\l/(r)|*(0)} wiU then also have 
a spectral representation which can be analysed 
with our methods. Fig. 2 is an example of such a 
procedure. 
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Figure 2. Spectral functions for 3 different super- 
positions of the ground-state and excited state 
smearings. 



This could be used to construct optimised 
smearing functions that contained Httle or no con- 
tamination from unwanted states. 

3.3. Error analysis 

An important issue is of course the error anal- 
ysis for MEM reconstructed spectral functions. 
One important point to note here is the fact that 
the spectral function is a density. It therefore 
does not make sense to assign errors to individ- 
ual points, since different points are correlated 
because e.g. the normalization of p is fixed. One 
can adopt an approach suggested in [pj in which 
the covariance of p around the maximum of Q 
is averaged over a. This is an approach within 
the logic of Bayesian statistics. We have taken 
a more reserved approach and asked how robust 
is the MEM prediction under a change of sam- 
ple. To investigate this we have created boot- 
strap samples from the original data and run the 
whole analysis for every such sample. For the 
ground-state smearing this procedure results in 
a stable prediction. The peak position does not 
change, only the shape of the peak. For the 1st 
excited state smearing the situation is still very 
good. The result of this analysis for the 2nd ex- 
cited state smearing is displayed in Figure 3. It 
is clear from this figure that one cannot trust the 
result for this smearing. 

Bootstrapped spectral function f (a))=p((o)/o)^ 




Figure 3. Spectral functions for 10 bootstrap 
samples for the 2nd excited state smearing show- 
ing the variation in shape and position of the 
peaks. The dotted lines are errors of the peak 
position based on 10 bootstraps. 



This shows that it is possible to investigate the 
statistical significance of MEM prediction by gen- 
eral bootstrap methods. 

4. Conclusions 

We have shown that with the Maximum En- 
tropy method one has access to the spectral func- 
tion of a current-current correlator which contains 
all the information about the excitation spec- 
trum of the theory in the channel represented by 
this current. We have also shown that one has 
good access to excited states when smeared oper- 
ators are used. In particular rotating correlators 
in a given basis of smearing functions, one can 
produce smearings which contain almost no con- 
tamination from unwanted states, complement- 
ing standard spectroscopy techniques. Since the 
method makes no assumptions about the shape 
of the spectral function, it can also be used at fi- 
nite temperature where little is known about its 
structure except for the absence of single particle 
(5-function peaks. We intend to use the method to 
investigate T and J /'^ melting above the critical 
temperature. 
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